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Abstract 



We work with two general factorization schemes in order to explore the conse- 
^ I quences of imposing scheme independence on g^{x,Q'^). We see that although the 

light quark sector is indifferent to the choice of a particular scheme, the extension 

>: 

0^ . of the calculations to the heavy quark sector indicates that a scheme like the MS 

^ ', is preferable. 

Q^ • The problem raised by the results of the EMC spin experiment [|1| was deeply influential 



on a substantial part of the 1990 's research in both theoretical and experimental hadron 



Oh. and particle physics. Their data implied that the quark singlet axial charge measured in 



a proton target, g^, was compatible with zero, while quark model calculations predicted 



k> ! fl'a to lie in the range of 0.6 - 0.7. After a series of experiments made at CERN, SLAG 

c^ I and HERA over the past 10 years, it is accepted today that g^ ~ 0.3, which is still far 

from those early theoretical expectations. 

In 1974, Ellis and Jaffe proposed a sum rule for the integral in x of gi{x), where 
they assumed that the sea quarks in the proton are not polarized. This implies that 
As, the helicity carried by the strange quarks, is zero. Experimentally, /q gi{x, Q^ = 
10 GeV^)dx = 0.120 ± 0.005 ± 0.006 ± 0.014 0, while the Ellis and Jaffe sum rule gives 
/o 9i{x, Q^ = 10 GeV^)dx = 0.176 ± 0.006. 

In the the parton model, the first moment of ^'^'(x) is given by Jq gi{x)dx = ^ga+^gl,+ 
^g^, with ga = Am — Ad the isotriplet axial charge, g^ = Au + Ad — 2As the octet axial 
charge, and g^ = Au + Ad + As. The parton model structure functions are, actually, QCD 
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structure functions with 0{a^) corrections. Beyond the parton model, the identification 
of the singlet axial charges with the sum of the quark helicities ceases to be true, because 
of the clash between a gauge invariant and a chiral symmetric renormalization procedure 
of the axial charge ^ ^ . 

The failure of the Ellis and Jaffe sum rule to agree with the experiments is translated 
into the non-equivalence between the quark singlet and octet axial charges. From the 
start there has been a large controversy on the mechanism responsible for g^ ^ Oa- ^^ 
the context of the parton model. As 7^ settles the question. However, as proposed by 
Altarelli and Ross [0, and Carlitz, Collins and Mueller |^, it is still possible to have As = 
if one takes into account the axial anomaly p] which appears in the QCD calculations of 
QilXjQ"^) at 0{as). Later, it became clear that these two different scenarios. As 7^ or 
an anomaly contribution, are simply related by a change of scheme defining the partons 
distributions and the coefficient functions. This will be, indeed, the main contribution of 
this work. We will argue that although the appearence of gluons in the first moment of 
gi{x, Q"^), in the light quark sector, is a matter of scheme preference, the introduction of 



heavy quarks suggests that a scheme where the gluons do not contribute, like the MS 
scheme, is preferable. 

A part of what is discussed in this work have already been adressed in the literature. 
Specifically, the importance in isolating the hard part of the photon-gluon cross section 
IP, |T0|, |ll|, |12[. However, some missconceptions still persist, mainly those connected with 
the heavy quark contribution to gi{x, Q"^), which is one of the motivations for the explicit 
discussions made here on the ways to calculate a polarized gluon coefficient function which 
is free of ambiguity in the infrared region. 

The choice of a factorization scheme is a refiection of the choice of a regulator and 
of a subtraction for the soft and collinear divergences appearing in the calculation of the 
partonic cross sections. In the specific case of the axial anomaly contribution to gf, much 
have been discussed about the ambiguity in the choice of a quark or of a gluon mass to 
regulate these divergences [|, |^ ^, |10[- In a satisfactory calculation, the hard part of 



the partonic cross sections should not present any ambiguity. The infrared singularities 



are present in the full virtual photon-gluon cross section |T^, |T^, Cg, and they appear 
explicitly when the Q^ — > oo limit is taken. As a general rule, the divergent (or soft) part 
of the cross section can be calculated from the expectation value of the quark singlet axial 
current between off-shell gluon lines ^, ||, |10|, p!5| : 



Aq\x,mlP',fi') = -2a,/i-(^-^) 



rf^-2fc, 



(2vr; 



D-2 



k'j_{l — 2x) — m 



D-2^ ' 



{kl + ml + x{l-x)P^y 



kl 



(/ci + m2 + x(l - x)P2)2 



;i) 



where P^ = —p^ is the gluon virtuality, m^ is the quark mass, and the number of quark 
flavors was set to 1. For an arbitrary number of flavors, Eq. (^ is multiplyed by Uf. 
The integral can be calculated in D dimensions as it stands, and the use of the modified 



Minimal Subtraction {MS) method to remove the UV divergence of Ag^ will define the 
coefficient functions and parton distributions in that scheme. A second option is to take 
from the start the limit Z) -^ 4, and use a cutoff /i^ to regularize the mass divergences. 
This is a momentum subtraction scheme, and the anomalous gluon contribution to the 
first moment of f^i (x, Q"^) will appeaifj. Explicitly, 



Aqf^x,mlP\^^') 



27r 



(2x - l)ln 



/i^ 



m„ 



^ml + x(l — x)P'^ ) ml + x(l — x)P2 



+ 1 



(2) 
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Ag^^(x,m,P^/x^) = ^ 



27r 

+ (l-x 



{2x - l)ln 



'/i^ + m^ + a;(l — x)P'^ 
m? + a;(l — x)P'^ 



fi^ 



2ml + a;(l - 2a;)p2 



(3) 



/i^ + m^ + x(l — x)P2 rrig + a;(l — x)P'^ 

Both Eqs. (^ and (^) are dependent on the rp?jP'^ ratio, which is not a real problem 

because they are only part of the gluon coefficient function. What configures a problem 
^Schemes like the AB of Ref. llfl], or JET of Ref. |0, belong to this class, and we will denote this 
class of schemes by [i schemes. 



is the attempt to draw conclusions about the possible anomalous gluon contribution to 
gi{x, Q"^) from those two equations. The standard procedure is to look at the subtracted 
partonic cross sections: 



C^'^x,Q\^^') 



Cg{x,mlP\Q^ 



Agf^(a:,m^,P^/i^) 



C^Ax,Q\^J?) = Cgix,mlP',Q') - Aq^ix,mlP',i^'). 



(4) 



The hard part of the cross sections should not depend on the m'^/P'^ ratio for Q'^,fi^ » 
m?, P^, which is the reason for the neglect of those two variables in the left hand side of 



Eq. (^. We also use the label MS to denote the fact that the usual MS gluon coefficient 
function is recovered in the large Q^ limit. The same for the Cj^ defined in a momentum 
subtraction scheme. 

In the region of low k±, the integrands of Eq. (|l|) and of Cg are equal. Hence, the fact 
that the RHS of Eqs. (^ turns out to be nonzero is a refiection of the large k± region, 
and of the UV regulator of Eq. (|l]): the redefinition of the parton distributions, through 
an absortion of finite parts of the cross section, is a matter of taste. Depending on the 
regulator chosen, one can also absorb or not the axial anomaly term into the redefinition 
of the parton distributions. Explicitly, when Q^,/i^ >> ■m'^,P'^, we have: 



cr(x,Q^/z2) 



as 
2^ 



(2x - i: 



MS-"(^ 



X 



- 1 



+ 2(1 -a;) 



«« 



C^{x,Q',^i') = -^(2x-l 
" In 



■4f)-"(^'- 



(5) 



Contrary to repeated claims in the literature |T^, |T^, [T^ , the schemes discussed here 
have a well defined separation of hard effects in the coefficient functions and soft effects in 
the parton distribuitions. In principle, the polarized light quark sector is well described 
by both of them. For a better understanding of both schemes, we show in Fig. || the 
integrals in x of Eqs. @ as a function of x, denoted by Iq^ = Jo C^'^^^{x, Q'^)dx. As is 
well known, /^ ^^^(x, Q^/x^) ^ g, and /q C^(x, Q^/i^) = -1, in units of a,/27r. The 



interesting feature is that the main contribution to both integrals conies from the large x 
region. In fact, x = 0.001 is already a good zero, while the x > 0.8 region is essential to 
give the integrals the value they have. 
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Figure 1: The integral from to x of the polarized gluon coefficient function as a function 



of X, for the MS and /x schemes, in units of as/2TT. schemes. 



The physical structure function is indifferent to which scheme is used to define the 
parton distribution and the coefficient functions. This is expressed as: 



J- n A/f c / N „q1\Tq , ^On J- 



+ ^g'','^{x)^C^'^{x,Q') + ^Ag^{x)^C^{x,Q'). (6) 



Although we did not write exphcitly the Q dependence of the various distributions, 
we remind the reader that only the singlet axial charge, in the MS scheme, has a Q^ 
dependent first moment. 

To 0{as), C,{x,Q^) = C^'{x,Q') = C^{x,Q') = 6{x - 1) + £^C«(x). Using 
Eqs. (|]), and the relation between the singlet axial charges between the two schemes, 
g^a'^ix, Q2) = gO,i^(^x, g2) _ ^^^n/(l - x) ® Ag{x, Q^), the second line of Eq. §) can 
be rewritten as: 



+ Ig'a^ix) ® Cl^ix, Q') + Ug^\x) ® Cf^ix, Q'), (7) 

where the term a;^/27r^(l — x) (X) Ag(x) Cj^^\x) was disregarded. We can now relate the 
remaining distributions and coefficient functions in the two schemes in the following way: 



C!;ix) = C'^^^'ix) + 6C,ix) 
Ag^ix) = Ag'^{x)+6g{x) 
^9li{x) + ^gl''{x) = ^9Y^{x) + ^9l^{x)+6q{x), (8) 

where 6Cg{x), 6g{x) and 6q{x) are some general functions, of 0{as)- Their specific form 
is not of interest to us at this given moment. However, the use of Eqs. (§) in Eq. (0), and 
the requirement that Eq. (Bl) is satisfied, produces the following consistency relations: 
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+6q{x) ® Cf^ix) + 5g{x) ® Cp(x) = 0, (9) 



+5q{x) ® C;(x) + 5(7(a;) ® ('c^^(a;) - ^r2/(l - a:)) = 0. (10) 

The first moment of Eq. @ certainly respects tlie equality, as /q Cg^^{x)dx = and 
/o Sq{x)dx = /o (5Cq(a;)(ia; = because of the conservation of the nonsinglet axial current 
^ The conservation of the nonsinglet axial current also imposes, from the first moment 
of Eq. (|10D, that Jq 5g{x)dx = 0, because Jq Cl^{x)dx = —nfas/2nAg. It follows that the 



first moment of the polarized gluon distribution is the same in the MS and /i schemes, 

independent of whether Ag contributes or not to the first moment of gi{x, Q"^), up to the 

{as/2'n-)'^Ag corrections we neglected before. This result is consistent with the fact that 

Ag starts contributing to gi{x, Q"^) at 0{as) only. 

Although both schemes are, in principle, equally good to describe f^i (x, Q^) in the light 

quark sector, we should also look at their behavior when heavy quarks are introduced. 

In particular, we do not want the hard part of the cross sections to depend on P^ once 

the mass of the heavy quark and the factorization scale are fixed. To investigate that, 

we calculate Eqs. @ as a function of P^ for the charm quark, with iric = 1.5 GeV 

and Q^ = fi^ = 10 GeV^. The resulting curves are shown in Fig. |^, normalized by the 

coefficient functions calculated with P^ = 0. It is clear that Cg^^ is independent of P^ in 

the range < P^ < ml. The same is not true for C'f, which shows a strong dependence 
onP2. 

Numerically, /^^ C^(a;, m^ = 2.25 GeV'^.Q'^ = 10 GeV'^)dx ^ 0.4, in units of 

Q^s(Q^)/2vr. Of course, this integral changes with Q'^, going to zero as Q^ -^ oo, but it is in- 
^ As the change of the coefficient function is dictated by the change of scheme of the anomalous dimen- 
sion, and the first moment of the nonsinglet anomalous dimension is zero due to current conservation, it 
follows that J„ 5Cq{x)dx = 0. 



dependent of P^ for fixed Q^. Hence there is a well defined contribution from gluons to the 



first moment oi gi{x, Q"^), in the MS scheme, which appears because of the relatively large 
mass of the charm quark. On the other hand, /g C^{x, m? = 2.25 GeV^, Q^ = 10 GeV'^)dx 
ranges from ~ —0.18, at P^ = 0, to ~ —0.135, at P^ = m^. Although the difference is 
not numerically significant ((0.18 — 0.135)(a;s/27r)A(7 ~ O.OOlA^f) as long as Ag is not 
very large, the use of the /i schemes is, to some degree, damaged. 

The inclusion of heavy quarks in the framework of perturbative calculation of structure 
functions have received great attention in the recent literature [0, ^, |^, |2^, |2^, ^^. 
These works have focused in the development of shemes that interpolate the pure photon- 
gluon fussion calculation from the region where Q^ ~ m|, to the usual massless approach 
(when Q^ >> ml). A fundamental point of these schemes is that the heavy quark 
is treated as a massless parton in the Altarelli-Parisi evolution equations, which will 
have Uf + 1 active flavours, while the quark mass dependence is fully kept in the graphs 
containing the heavy quark lines in the calculation of the coefficient functions. These 
schemes are generally refered to as Interpolating Schemes (IS). 

The coefficient functions in Eq. (^ incorporate the full mass corrections, and are 
reduced to the massless case in the limit of large Q^. Hence, they are suitable for the 
calculation of the polarized structure functions for Q^ ~ m| and Q^ » ml, in the spirit 
of the IS. As in the IS the light and the newly introducded heavy quark distributions 
should be defined in the same scheme, and as the calculation of the C^ for a heavy quark 



is ambiguous, it follows that, strictly, the MS scheme is formally superior to the /i scheme 
for the calculation of gi{x, Q"^). 

As a last remark, we want to stress that the amount of polarized heavy quark in the pro- 
ton is not given by the integral in x of Eqs. @) and (^, or from the integral of Cg(x, Q^) for 
a given quark mass. From them, one would conclude that /q Aq^^^{x,m'l,P'^,jj?)dx = 
for /i^ = 771^ >> P^, while /q Ag^(x, m^, P^,yU^)(ia; = Q;s/47r in that same limit. In a 
framework where heavy quark mass effects are systematically included, one should intro- 



duceQ a polarized heavy quark distribution in the proton (A/i(x, Q^)) at the factorization 



scale /i^, with A/i(x, /i^) = 00. As we saw here, both MS and // schemes are suitable 



for this purpose once Eqs. (H) are given, although, in principle, the MS scheme has the 
technical advantage of having a P^ independent gluon coefficient function in the heavy 
quark sector. 

I would like to thank X. Ji and A. W. Thomas for valuable discussions. This work was 
supported by FAPESP (under contracts 96/7756-6 and 98/2249-4). 

References 

[1] EMC Collaboration, J. Ashman et al, Phys. Lett. B 206 (1988) 364; Nucl. Phys. 
B328 (1989) 1. 

[2] J. Ellis and R. L. Jaffe, Phys. Rev. D9 (1974) 1444. 

[3] SMC Collaboration, B. Adeva et al, Phys. Rev. D 58 (1998) 112002. 

[4] W. Bardeen, Nucl. Phys. B 75 (1974) 246; R. J. Crewther, Acta Physica Austriaca 
Suppl. XIX (1978) 47. 

[5] G. Veneziano, Mod. Phys. Lett. A4 (1989) 1605; S. D. Bass, R. J. Crewther, F. M. 
Steffens and A. W. Thomas, Piep-ph/9701213| . 



[6] G. AltareUi and G. G. Ross, Phys. Lett. B212 (1988) 391. 

[7] R. D. Carlitz, J. C. Collins and A. H. Mueller, Phys. Lett. B214 (1988) 229. 

[8] S. L. Adler, Phys. Rev. 177 (1969) 2426; J. S. Bell and R. Jackiw, Nuovo Cim. 60 
(1969) 47. 



''if the heavy quark contribution to gf (.t, Q^) is calculated through the photon-gluon fussion proccess 

only, and its higher order corrections, a polarized heavy quark distribution is never introduced. 

^It is assumed that there is no intrinsic heavy quark polarization in the proton. See Ref. 123] for a 

different point of view. 



[9] G. T. Bodwin and J. Qui, Phys. Rev. D41 (1990) 2756. 

[10] L. Mankiewicz, Phys. Rev. D43 (1991) 64. 

[11] F. M. Steffens and A. W. Thomas, Phys. Rev. D 53 (1996) 1191. 

[12] Hai-Yang Cheng, Int. J. Mod. Phys. All (1996) 5109. 

[13] S. D. Bass, N. N. Nikolaev and A. W. Thomas, Report No. ADP-133-T-80, 1990 
(unpubhshed) 

[14] W. Vogelsang, Z. Phys. C50 (1991) 275. 

[15] G. P. Lepage and S. J. Brodsky, Phys. Rev. D22 (1980) 2157. 

[16] R. D. Ball, S. Forte and G. Ridolfi, Phys. Lett. B378 (1996) 255. 

[17] E. Leader, A. V. Sidorov and D. B. Stamenov, Phys. Lett. B445 (1998) 232. 

[18] E. Leader, J. Phys. G 25 (1999) 1557. 

[19] M.A.G. Aivazis, J.C. Collins, F.I. Olness and W.-K. Tung, Phys. Rev. D 50 (1994) 
3102. 

[20] F.M. Steffens, Nucl. Phys. B 523 (1998) 487. 

[21] S. Kretzer and I. Schienbein, Phys. Rev. D 58 (1998) 094035. 

[22] M. Buza, Y. Matiounine, J. Smith and W.L. van Neerven, Eur. Phys. J. C 1 (1998) 
301. 

[23] R. S. Thorne and R. G. Roberts, Phys. Rev. D 57 (1998) 6871; Phys. Lett. B 421 
(1998) 303. 

[24] J. C. Collins, Phys. Rev. D 58 (1998) 094002. 

[25] S. D. Bass, S. J. Brodsky and I. Schmidt, |hep-ph/9901244 . 



10 



1.5 



1.25 



o 
II 

CVl 

CL 

OJ 

O 

o 

CM 

CL 

CM 

o 
o 



0.75 



0.5 



)j. scheme 

Mod. MIn. Sub. Scheme 



0.5 



1.5 



P'(GeV') 



Figure 2: The integrated (in x) polarized gluon coefficient function in the MS and /i 
schemes, as a function of P^. 
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